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The Coulomb problem for vector bosons W ± incorporates a well known difficulty; the charge of 
the boson localized in a close vicinity of the attractive Coulomb center proves be infinite. This fact 
contradicts the renormalizability of the Standard Model, which presumes that at small distances 
all physical quantities are well defined. The paradox is shown to be resolved by the QED vacuum 
polarization, which brings in a strong effective repulsion that eradicates the infinite charge of the 
boson on the Coulomb center. This property allows to define the Coulomb problem for vector bosons 
properly, making it consistent with the Standard Model. 
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I. INTRODUCTION 

Consider a charged vector boson, which propagates in 
the Coulomb field created by a heavy point-like charge 
Z assuming that the boson is massive, its mass being 
produced via the Higgs mechanism; the IT ± -bosons give 
an example. We study relativistic effects in this Coulomb 
problem. A situation where they can be important arises, 
for example, for small primordial charged black holes 
since an impact of their Coulomb field on a W-boson 
prevails over the gravitational field. 

It has “always” been known that there is a difficulty 
in the Coulomb problem for vector bosons. Soon after 
Proca formulated theory for vector particles [j] it be¬ 
came clear that it produces inadequate results for the 
Coulomb problem HIS This fact inspired Corben 
and Schwinger [B] to modify the Proca theory, tuning 
the Lagrangian and equations of motion in such a way as 
to force the hyromagnetic ratio of the vector boson to ac¬ 
quire a favorable value g = 2. Later on the formalism of 
[5] was found to have a connection with the non-Abelian 
gauge theory Q, which makes it relevant for the present 
day studies. A role of the identity g = 2 was thoroughly 
discussed in literature, see e. g. Ref.jJ .8]. 

Ref. H found a realistic discrete energy spectrum for 
the Coulomb problem for vector bosons. However, it dis¬ 
covered also a fundamental flaw in the problem. For two 
series of quantum states the charge of the vector boson 
located on the Coulomb center turns infinite, which indi¬ 
cates the fall of the boson on the center. One of these se¬ 
ries has the total angular momentum zero, j = 0, another 
one has j = 1 (being further specified by a label “7—3/2”, 
see Section nvn . This effect takes place for arbitrary 
small value of the Coulomb charge Z, which is physi¬ 
cally unacceptable. Moreover, it takes place at small dis- 
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tances, while the renormalizability of the Standard Model 
Ref-0 guarantees that there should be no problems of 
this type. All this indicates that the Coulomb problem 
is poorly defined. Moreover, there exists a contradiction; 
the Coulomb problem derived from the Standard Model 
produces results, which challenge the Model itself. 

This difficulty was inspirational for several lines of re¬ 
search. Early efforts are summarized in Ref. E3- More 
recent Refs. n, El E3 suggested a new, refined mod¬ 
ification of the formalism for vector bosons. Ref. 0 
claimed that it complies with results of Corben and 
Schwinger. Some authors considered other forms of the 
equation governing vector bosons 0E1E3, which pro¬ 
duce more acceptable results for the Coulomb problem, 
but this advantage is partially undermined by the fact 
that it does not step from a renormalizable theory. 

However, in spite of a progress made over the years, 
there still exists a contradiction between the difficulty in 
the Coulomb problem for vector bosons and the renor¬ 
malizability of the Standard Model. We find a clear way 
to resolve this contradiction, formulating the Coulomb 
problem for vector particles properly, within the frames 
of the Standard Model. Our main observation is that 
the polarization of the QED vacuum has a profound im¬ 
pact in the problem forcing the density of charge of a 
vector boson to decrease at the origin, thus making the 
Coulomb problem stable, well defined. This decrease has 
an exponential character for the j = 0 state. For the 
j = 1,“7 — 3/2” state the suppression is of a power-type. 
In both these states the suppression eradicates the diffi¬ 
culty of the Coulomb problem. 

From the first glance this result looks surprising. Pre¬ 
sumably, the vacuum polarization is meant to make the 
attractive Coulomb field only stronger, which should re¬ 
sult in an increase of the charge density at the origin. In 
addition to this, the vacuum polarization for spinor and 
scalar particles in the Coulomb field is known to produce 
only small, perturbative effects. In contrast, we claim a 
strong reduction of the charge density for the vector par¬ 
ticle. To grasp a physical mechanism involved it is nec- 
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essary to notice that the equation of motion for vector 
particles incorporates a particular term, which explicitly 
depends on the external current and has no counterparts 
for scalars and spinors, see the last term in Eq. OB . 
Precisely this term brings in a strong effective repulsion, 
which stems from the vacuum polarization and makes the 
Coulomb problem stable, well defined. 

The renormalizability of the Standard Model means 
that if all essential processes are taken care of, then the 
infinite charge of a vector boson located at the Coulomb 
center is eliminated. It is known that the amplitude of 
the photon exchange between leptons or/and quarks at 
high transferred momenta should be considered alongside 
exchange by the Higgs and -Z-bosons. From this perspec¬ 
tive the catastrophic behavior of the charge density of a 
vector boson at small distance, i. e. at large transferred 
momenta, in the Coulomb problem could have been con¬ 
sidered as an indication that the Coulomb problem for 
vector bosons should include the processes related to the 
Higgs and .Z-bosons exchange from the very beginning. 
In contrast to this widely spread presumption we find a 
way to formulate the Coulomb problem for vector bosons 
entirely in terms of the W and electromagnetic fields, as 
a pure QED problem. 

A complete Standard Model calculation, where all pos¬ 
sible processes are accounted for accurately, would re¬ 
quire specific information on the nature of a heavy par¬ 
ticle that creates the Coulomb field, i. e. on all its quan¬ 
tum numbers related to the Standard Model. This infor¬ 
mation is not necessarily feasible. A simple example give 
primordial black holes; it is not easy to assert with cer¬ 
tainty whether they have, or have not the weak charge, 
and what are their other quantum numbers in the Stan¬ 
dard Model. Same questions arise in relation to other 
possible candidates for the heavy Coulomb center. As a 
result, a presumption that the exchange of the Higgs and 
.Z-boson should play a basic role in the Coulomb problem 
leads to complications. It is fortunate therefore that the 
detailed information on properties of the heavy particle 
proves be redundant, that the Coulomb problem can be 
properly defined using the only physical parameter of a 
heavy particle, its electric charge. 

This point of view, which is advocated in the present 
work keeps the Coulomb problem simple and transpar¬ 
ent. On the other hand, it also allows one to include all 
other processes, which are left outside the scope of the 
Coulomb problem, by means of perturbation theory. Our 
preliminary calculations indicate that the exchange by 
the Higgs and .Z-bosons, as well as possible processes with 
lepton or quark exchange, give only small corrections. 
The reason stems from the fact that the found wave func¬ 
tions for vector bosons are suppressed at small distances. 
Consequently, the small-distance processes with the ex¬ 
change by Higgs and .Z-bosons are also suppressed (the 
exchange by a lepton or quark contains the vanishing at 
the Coulomb center fields, which describe the IT-boson). 

In section QH the Corben-Schwinger formalism for 
charged vector bosons is derived directly from the Stan¬ 


dard Model. The pure Coulomb problem is discussed 
in Sections HillVI and several Appendixes. This analy¬ 
ses follows Ref. j3, but some important details, including 
the non-relativistic limit (Section IlIIll and the eigenvalue 
problem for j = 0 states (Section 0 are discussed in 
more detail. Sections I V II V 111 present the main result of 
the paper. They show that the QED vacuum polarization 
plays a defining role in the problem, as was first noticed 
in our previous work 0- The units h = c = 1, e 2 = 47ra 
where e < 0, are used below. 

II. W-MESONS IN ELECTROMAGNETIC 
FIELD 

A. IT-bosons in Standard Model 

Consider boson fields in the electroweak part of the 
Lagrangian of the Standard Model, see e.g. Ref. fl9j. 

£ = ~~ ^ (d^Av — d^A^ + g A^ x A v ) , (2-1) 

- d v B^f + i D^+D^ . 

Here A^ and B M are the triplet of SU( 2) and the U{ 1) 
gauge potentials respectively (abridged notation is used 
here). The covariant derivative takes into account 
that the Higgs field $ has a hypercharge Y = 2, which 
describes its interaction with the 17(1) field, and is trans¬ 
formed as a doublet under the SU (2) gauge transforma¬ 
tions. Taking the unitary gauge one can present it via 
one real component 

* = ( J ) , <t> = <t> + ■ (2-2) 

Assuming that the scalar field develops the vacuum ex¬ 
pectation value <f> = (j>o and the Higgs mechanism takes 
place, one finds that the gauge held can be presented as a 
new 17(1) held A M , and a triplet of massive helds W±, Z^ 

A M = — sin 9 + cos 9 B^ , (2.3) 

Z^ = cos 9 A^ + sin 8 B^ , (2.4) 

W, = (Al-iAl)/V2, (2.5) 

Here = W~ represents the VF-boson with charge 
e = — |e|, and 6 is the Weinberg angle. 

Expanding the Lagrangian Eq. JO in the vicinity of 
</> = </>o and retaining only bilinear in the helds W \ L , W+ 
terms, including their interaction with the electromag¬ 
netic held, one derives an effective Lagrangian 

C w = -i(V„W„ - V 1/ W , M ) + (V'W 1 '- V U W^) 

+ie W+W v + to 2 W+W^ , (2.6) 

which describes the propagation of W-bosons in an exter¬ 
nal electromagnetic held. Here m is the mass of W. The 
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external field is accounted for in Eq. (El) in the deriva¬ 
tive + ieA „ and by the term with the field 

F^ — d v A^. The first and the last terms in 

Eg. (12. (ID are present in the Proca formalism T|, while 
the second one was introduced by Corben and Schwinger 
@- 

From Eq. ESI) one derives the classical Lagrange equa¬ 
tion of motion for vector bosons 


Using this result, the current Eq. EH} can be written in 
a compact form 

jjf = -ie(w+V M fU" + 2V !/ W+VP v -c.c.) 

-fr {W+W v + w+w^f . (2.16) 

771 z p 

Here c.c. refers to two complex conjugated terms. 


(V 2 + to 2 ) W" + 2 ie W v - W v = 0 . (2.7) 

Here an identity [V^V,,] = ieF was used. Taking a 
covariant derivative in Eq. E3 one finds 


m 2 V /1 W’ ,i + iej^W* = 0 , (2.8) 


where 


f = d v F v » . (2.9) 

is the external current, which creates the external field 
F l 'U Evaluating from Eo. ll2.8D and substituting 

the result back into Eq. E3 one rewrites the latter one 
in a more transparent form 

(V 2 + to 2 ) W* + 2 ieF^Wu + ^W{j v W v ) = 0. (2.10) 

m z 

This equation of motion for vector bosons was suggested 
in Ref. [a]- The coefficient 2 in front of the second term 
ensures that the g-factor of the boson takes the value 
g = 2, see Eq. Em below. 

The derivation outlined shows that Eq. EE3> repre¬ 
sents the classical equation of motion for W-bosons in 
the external electromagnetic field, which is valid within 
the frames of the Standard Model. This equation has 
similarities with the Klein-Gordon and Dirac equations 
(if the latter one is written as the second-order differen¬ 
tial equation), but there is also an important distinction. 
It is produced by the last term in Eq. Em which ex¬ 
plicitly contains the external current; there is no similar 
terms for scalars and spinors. We will see how important 
this term is, when we discuss the vacuum polarization. 

We will use below a current of vector bosons j ^, which 
can be obtained by considering a variation of the La- 
grangian Eq. ES> under variation of A which yields 


W 
J11 

= if+if+if, 

(2.11) 

j (1) 

J t* 

= -ie(W+V ll W v -V li W+W v ) 

(2.12) 

J (2) 

J fl 

= -ieiVW+W 1 ' 

(2.13) 

7 (3) 

J n 

= -ieriW+Wv-W+Wn) . 

(2.14) 


Differentiating in Eq. eh term by term and taking into 
account Eq. El one verifies that 

if = jf ) -ie(W+V l/ W v -V , 'W+W fi ) (2.15) 

= if - W Z W " + } ^ ■ 


B. Static electric field 


Consider a static electric field described by the electric 
potential A 0 = A 0 (r) and charge density p = p(r) = 
— AA 0 . For a stationary state of the W-boson one can 
presume that 

V 2 W„ =-(e - UfW^ , (2.17) 

where e is the energy of the stationary state, and U = 
U(r) = eAo is the potential energy of the TU-boson in 
the electric field. Eq. El in this case gives 

w = (e-U - T) -1 V • W. (2.18) 


The four-vector = (Wo,W) is presented here via 

the three-vector W and the modifies zeroth-component 
w = iWq. In order to simplify notation we introduce also 
a very important for us quantity T = Y(r), 


ep A U 

m? to 2 


(2.19) 


Eqs. (12.8D , Ii2.17[) show that this definition complies with 
EH - The quantity T appears in the equations of mo¬ 
tion alongside the initial potential C7, see e.g. Eq. EH - 
In this sense it plays a role of an effective potential en¬ 
ergy, which is specific for vector bosons. We will call it 
the Y-term, or Y-potential. In this notation Eq. EH 
reads 


((e - Uf- m 2 )W = -AW - 2Vt/w - V(Yw), (2.20) 
((e — U) 2 — to 2 )w = —Aw + 2V/7 • W (2.21) 

+{e-U) Yw. 

A relation between w and W given by Eq. eh shows 
that among four equations of motion Eqs. EH, EH 
only three are independent, precisely what one expects 
for massive vector particles. 

It will be useful to present Eq. ESI in a slightly dif¬ 
ferent form, which can be derived by combining it with 
Eci. (l2.18D and using an identity AW = V X (V X W) — 
V(V- W), which gives 

((e — U) 2 — m 2 )W = V X (V X W) (2.22) 
— (e — t/)Vw — VI/ w . 

From the expression for the current of vector bosons 
Eq. EH one derives the charge density 

p w = 2e^(e — U)(W + - W + w + w) 

+ W + -Vw + W-V w + — Y w + w 


(2.23) 
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C. G-factor 


$ of the vector boson 


The behavior of vector bosons in the homogeneous 
magnetic fields was studied in detail, see e.g. [lCj and 
references therein. The spectrum of this problem reads, 
see Section El 

£ 2 = m 2 + pi + 2\e\B (n + 1/2 + a) . (2.24) 

Here n = 0,1... specifies the Landau levels, and a = 
— 1,0,1 gives a projection of spin S' = 1 of the vector 
boson. Ea. 12.2411 shows that vector bosons possess the 
magnetic moment 


p, = eS/m , (2.25) 

which means that the magnetic ^-factor is <7 = 2. 


III. NON-RELATIVISTIC LIMIT 

Consider a vector boson in a static electric field with 
the potential energy U = eA^(r). If we presume that 
the non-relativistic approach is valid, which needs that 
\U\ -C to, then in the lowest order of the perturbation 
theory in powers of U/m one immediately finds from 
Eos. 12.2011 . 12.1811 


EW = ——A W + UW. (3.1) 

2 TO 

Here E ~ e — to is the energy, the vector W plays a 
role of the wave function for the vector boson, and the 
non-relativistic Hamiltonian on the right-hand side has a 
usual form for a massive charged particle. 

Let us find corrections to Eq. El induced by rela¬ 
tivistic effects. The wave function of the massive vector 
particle $ is well defined in the rest frame. Therefore 
the vector W, which describes the moving vector par¬ 
ticle, inevitably deviates from the wave function A 
relation between W and <1? is easy to articulate for the 
free motion, when it is given by the Lorentz boost, see 
e.g. a book f20| , 


W = $ + 


P(P $) 

m(m + e) 


(3.2) 


Generically, the potential energy brings in complications, 
but within the necessary accuracy we can neglect them, 
presuming also that e ~ to. Then Eq. 13.211 gives 


W ~ $ + 


P(P&) 

2 TO 2 


(3.3) 


where p = — iV. This relation plays a role similar to the 
Foldy-Wouthuysen transformation |2lJ for fermions. 

Substituting Eq. 13.31) in Eas. 12. 21)1) . 12.181) and expand¬ 
ing the latter ones in powers of U/m one finds the fol¬ 
lowing Schrodinger-type equation for the wave function 


Ha 


SHij 


(iL + U ) 6ij + 6Hij ’ 

P A , _ F ■ (p x Sjj) 
8to 3 u 2 m 2 

6 to 2 dridrj Sij'j 


A U 
6to 2 u 


(3.4) 

(3.5) 

(3.6) 


Here i. j = 1,2,3 label components of three-vectors, S 
is the spin, which operates on a vector V according to 

E'ij Vj = fejjfcVfc. 

The relativistic correction to the Hamiltonian 5H of 
vector particles is given in Eq. ESI- It is instructive to 
compare this correction with the known Darwin Hamilto¬ 
nian SHd , which accounts for relativistic effects for spinor 
particles 


SH d 


p 4 F ■ (p x s) A U 
8 to 3 2 to 2 8 TO 2 


(3.7) 


Here s = cr /2 is the operator of spin for spinor parti¬ 
cles. The three terms in the first line of Eq. 13.61) resemble 
their counterparts in Eq. (1X71) . the only distinction is the 
numerical coefficient in front of the term with A U. We 
conclude that these three terms have conventional mean¬ 
ing, describing the relativistic correction to the kinetic 
energy, the spin-orbit interaction, and the contact cor¬ 
rection to the potential. The coefficient in front of the 
term responsible for the spin-orbit interaction in Eq. 13.51) 
complies with the hyromagnetic ratio g = 2 of the vec¬ 
tor boson, if one presumes that the Thomas “one-half 
rule” is applicable for vector particles the same way as 
for spinors. 

The last, forth term in Eq. E3> finds no counterpart 
in the Darwin Hamiltonian. It is instructive to write a 
contribution of this term to the energy shift 

I I" r ) 2 A 

SE ° = < 38 > 

Qij = • (3-9) 

Eas. 13.81 ). 13.91) show that plays a role of the density of 

the quadrupole moment for vector bosons. We conclude 
that the last, forth term in Eq. 13.61) indicates that vector 
bosons have a quadruple moment. 

From the first glance the contact and the quadrupole 
terms in the Eq. EH have similarity with the T-term 
in Eq. EH- However this resemblance is coincidental, 
since the T-term does not contribute to ESI, which takes 
into account corrections of the order of ( Za ) 2 . En. 12.1 9|) 
allows to estimate the Y-potential as T ~ ( mro)~ 2 U ~ 
(. Za) 2 U , where ro = ( Zam ) _1 is the Bohr radius. The 
T-potential comes into the equation of motion with the 
factor w, see the last term in Eq. a - Eq.Eni gives 
an estimate w ~ (mro)~ 1 \W\ ~ Za\W\. Overall, an 
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estimate for the correction produced by the Y-term in 
Eg. (12.2011 is ~ ( Za ) 3 , which means that the T-terrn is too 
small to contribute to Eq. iTHl . Thus, the contact and 
quadrupole interactions in Eq. (13.01) have no direst con¬ 
nection with the Y-term. This fact makes a difference in 
coefficients in front of the contact term in Eq. (EZ§ and 
the Y-term in Eq. (mm acceptable. In particular, the 
fact that they have opposite signs produces no contra¬ 
diction. 


IV. COULOMB PROBLEM 


vectors, , Y}™'* defined conventionally, see 

and Appendix [C] Generically, one can present the vec¬ 
tor W as a linear combination of three spherical vectors 
with the given value of j. It is convenient to refer to the 
three terms in this combination as the electric, longitu¬ 
dinal and magnetic modes (or polarizations) of a vector 
boson. The parity conservation simplifies the problem 
further on. The state with the magnetic polarization, 
which parity is different from the parity of other two 
modes, is not coupled with these modes. Therefore the 
magnetically polarized mode can be written in a simple 
form 



Consider the pure Coulomb held, presuming that it is 
created by a point-like heavy object with charge Z > 0. 
Then for r > 0 one has 

[/ = -— , T = 0 . (4.1) 

r 

The second identity here follows from Eq. m> . 


A. Perturbation theory 


Let us treat the Coulomb problem using the non- 
relativistic perturbation theory. Take the non-relativistic 
Ea. 113 .111 as a starting point, and consider the Hamilto¬ 
nian Eq. 13. (ill as a perturbation. Conventional calcula¬ 
tions, see Appendix m lead to the following result for 
the shift of the energy level characterized by the main 
quantum number ro, orbital momentum l and total angu¬ 
lar momentum j = l, H 1 


fiEnlj — 


m(Za) 4 




(4.2) 


This formula is similar to the one that describes the en¬ 
ergy shifts for spinor particles; the only distinction comes 
from values of j in Eq. (TO) . which are integers for vec¬ 
tor particles and half-integers for spinors. The order of 
several lowest levels shows the following pattern 


W = fY (4.4) 

where f = f(r) is the radial function. The two modes re¬ 
lated to electric and longitudinal polarizations have same 
parity, which makes coupling between these modes pos¬ 
sible. One needs therefore to consider them on the same 
footing assuming that 

W = uY$+vY® , (4.5) 

where u = u(r), v = v(r). We will refer to them as 
electro-longitudinal modes, or polarizations. 


C. Magnetic polarization, j > 1 

For the magnetic mode the angular momentum is re¬ 
stricted j > 1 (the magnetic spherical vector is not de¬ 
fined for j = 0, see Eq. CD)- Substituting Eq. (IQI) into 
Eq. CZP one finds the following equation for the radial 
function / 


(A j + (e + Za/r) 2 — m 2 ) / = 0 . 


Here A j is 


A? = 


1 d 
r 2 dr 



j{j + 1 ) 

r 2 


(4.6) 


(4.7) 


n = 1 lsi; 
n = 2 2p 0 , 

{2a!, 

2pi}, 

2p 2 ; 

The form of Ea. (14.61) coincides with the Klein-Gordon 
equation. Therefore the spectrum of the magnetic mode 
' ' ' replicates the spectrum of scalar particles, which is given 

O 

& 

CO 

II 

fi 

{3si, 

3pi, 

3di}, 



3^2}, 3d 3 . 

by the Sommerfeld formula 


Here the atomic-like notation nlj is adopted, the brackets 
combine together the degenerate energy levels. 


e = m 



(Za) 2 \ 

(l + n-j- 1/2) 2 J 


- 1/2 


(4.8) 


B. Central field 

Consider the static central electric field (the Coulomb 
problem gives an important example). The conservation 
of the total angular momentum j in this field allows one 
to separate the angular variables. We will use for this 
purpose the electric, longitudinal and magnetic spherical 


Here 

7= ((j + l/2) 2 -(Za) 2 ) 1/2 . (4.9) 

In Ea. 14.81) n = 1,2 ... plays a role of the main quantum 
number. In the non-relativistic limit the magnetic mode 
corresponds to the states 2p\, 3^2,4/ 3 ,.... 
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D. Electro-longitudinal polarizations, 3 > 1 


Consider electro-longitudinal polarizations, when the 
vector W is given by Eq. 14.511 . Substituting it into 
Eos. 12.2011 , 12.1811 and using the properties of the spherical 
vectors from Appendix El one finds a system of coupled 
equations for radial functions u, v 

(A j + (s + Za/r) 2 - m 2 )u = -2 y/j(j + 1) , (4.10) 

(Aj + (e + Za/r ) 2 - m 2 )v = -2 y/j(j + 1) (4.11) 


1 r 2 j.2 

Here w = w(r) denotes the radial part of w. Using 
Ea. (l 2 ~T 8 l) one finds for it 


w = w Y jm , 


w = 


1 

e + Za/r 



(4.12) 

.(4.13) 


Eqs. CUES-ED are sufficient to define the functions 
u, v, but it is convenient to compliment them by the ra¬ 
dial form of Eq. (irm which reads 


(A j + (e + Za/r) 2 — m 2 )w = —-— . (4-14) 

Let us verify first that Eqs. cm cm describe two dif¬ 
ferent modes. Consider with this purpose distances so 
small that m <C Za/r , where the potential energy dom¬ 
inates over mass. In this region Eqs. cm cm reduce 
to 


dr 2 


dr 2 


2 d 
r dr 


4 d 
r dr 


(Za) 2 -j(j + l)\ 

-- ) U 

(Za) 2 -j(j + l) + 2\ 
r 2 J 


(4.15) 



v = 0. (4.16) 


One derives from Eqs. 14.15H . 14.Hill that there exists a 
mode, in which at small distances v is small, |u| <C |u|, 
which means that in this region the polarization is pre¬ 
dominantly electric. From Eq. n~T5i) one finds that this 
mode satisfies the following asymptotic conditions at 
r —> 0 


u->ar 7 172 , |u| -C |it| . (4-17) 

We will call it the “7 — 1/2” mode below. 

In order to find the second mode let us assume the 
following asymptotic behavior for r —> 0 

u->br v , (4.18) 

v —> cr 1/ , (4-19) 

Substituting Eas. 14.1811 . 14.191) in Eqs. 14.15(1 , 14.Kill one 
finds a system of two homogeneous linear equations, in 


which v plays a role of the eigenvalue. Solving this sys¬ 
tem one finds v and the ratio c/ 6 , deriving 

u -> 6 r 7-3/2 , (4.20) 

v -> 6 7 ~ 1//2 r 7-3 / 2 . (4.21) 

VJU+T) 

This mode will be referred to as the “7 — 3/2” mode j22^ . 

Let us find now the discrete energy spectrum. Intro¬ 
duce a function g = g(r) 


g = Zau + y/j(j + l)i 

= Zau+ VjU r7 + } ) 
£ + Za/r 


/ . / . TT 4T (I'D 

~VdU + 1)7 + -/// + 


(4.22) 
2v' 


Here Eq. IXT5I ) was used in the second identity. Taking 
the corresponding linear combination of Eqs. 633 , 631 
one finds that g satisfies the Klein-Gordon equation 


(Aj + (e+ Za/r) 2 — m 2 )g = 0 . (4-23) 


This result leaves only two options; either g equals zero 
identically, or, alternatively, the spectrum of electro¬ 
longitudinal modes can be found from Eq. 14.231) . The 
first alternative takes place for j = 0 , when only the lon¬ 
gitudinal mode is present. The function u in this case 
should be taken as zero, which makes zero also the func¬ 
tion g in Eq. 14.2211 . Thus, Eo. 14.231) provides no help for 
j = 0 states. 

For j > 1 the function g is nonzero, for both “ 7 — 1/2” 
and “7 — 3/2” modes, see Appendix IU1 Eo. 14.231) defines 
the spectrum, which therefore satisfies the Sommerfeld 
formula Eq. In the non-relativistic limit the mixed 

electric-longitudinal modes correspond to the following 
states with j > 1: lsi, 2p2, 3di, 3^3, 4/2 .... 


E. Longitudinal polarization, j = 0 


Consider zero angular momentum j = 0, which corre¬ 
sponds to purely longitudinal polarization, see Eq. ED- 
The state with j = 0 is described by one radial function 
v = v(r), 


W = v n , n = r/r . (4.24) 

The radial function v satisfies Eqs. in which 

the function u is to be put to zero (electric polarization 
for j = 0 is impossible). These equations therefore yield 

d 2 v 2 dv (, „ . . 2 

■ 3-7 -t-r- + \{e + Za/r)" - m 2 ) v (4-25) 

dr * r dr \ / 

2v 2 Za 1 / dv 2?A 

r 2 r 2 e + Za/r \dr r J 


In order to make the physical meaning of this equation 
more transparent let us eliminate the first derivative by 
means of a substitution v —> tp 


Za 


£ + 


Za 


r 


1 


■V 


(4.26) 
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where it is convenient also to scale the radial variable 
r —> x 


Za 

r =- x , 

£ 


(4.27) 


assuming ip = p(x). In this notation Ea. 14.2511 can be 
rewritten as a conventional Schrodinger-type eigenvalue 
problem 


H p = —x 2 p , 

d 2 2 {Za) 2 (Za) 2 

dx 2 x x 2 


2 

(a: + l) 2 


(4.28) 

(4.29) 


Eqs. run .rum give then the function p, while 
transform it into v and W. The function 
p exhibits the following behavior at the boundaries 

tp oc exp(— xx) , x —> oo , (4.38) 

p oc x 1+1 / 2 , x —> 0 . (4.39) 

Ea. 14.3611 gives the spectrum 

( \ _1 / 2 
1 H-----] , n = 2, 3 ... (4.40) 

(7 + n — 1/2) 2 J 


where — >r 2 , which plays a role of an eigenvalue, is related 
to the energy of the discrete level 

= (Za) ^ > 0 . (4.30) 

The operator H in Eq. lTP5t possesses three singular 
points, x = 0, x = 00 and x = — 1. The last one lies 
in the non-physical region, but it presents an obstacle 
for an analytical study anyway. One can overcome this 
difficulty using a substitution tp —> p 

<P = (4~ + (7 + I/ 2 ) -Tt) ^ ' ( 4 - 31 ) 

\dx x x + 1 J 


It can be shown that tp satisfies an eigenvalue problem 


Hip = —x 2 p , 

d 2 (Za) 2 

dx 2 x 


(4.32) 


+ (i+))(t + |)^. (4.33) 


The main result of the transformation Eq. (TOTli is that 
the operator H has only two singular points, x = 0 
and x = 00 . An interesting method, which allows one 
to “invent” the substitution Eq. (TOD) and derive then 
Ea. 14.3211 is presented in Appendix lEl It takes its origins 
in an elegant treatment of quantum mechanics developed 
by the Gotingen School and known as matrix mechanics. 

A regular at r = 0 solution of the eigenvalue problem 
(TO!! reads 


p = e~* x x L+1 F ( L + 1 - 


(Za) 2 


,2L + 2, 2xx . (4.34) 


Here F(a,(3,z) is the confluent hypergeometric function 
and L is defined by 


L = 7+1/2. (4.35) 


To make the solution given by Eq. (TOl regular at infin¬ 
ity one should assume that 


K = 


(Za) 2 
L + n — 1 


(Za ) 2 

7 + n — 1/2 ’ 


n = 2,3... , (4.36) 


The corresponding eigenfunctions are given by Eq. 623J), 
in which the hypergeometric function is reduced to a 
polynomial 

tp = e~ ><x x' 1+3/2 F(2 -n, 2 7 + 3, 2xx) . (4.37) 


which complies with the Sommerfeld formula Eq. 14.811 . 
In the non-relativistic limit the longitudinal mode corre¬ 
sponds to the following states with j = 0 : 2 p 0 , 3po> 4po ■ ■ • 


F. Summary for Coulomb problem 

Our discussion of the Coulomb problem for vector par¬ 
ticles confirms that for all polarizations and all angular 
momenta j the discrete energy spectrum is described by 
the Sommerfeld formula Eq. i TOI i. as was first found by 
Corben and Schwinger §5]. 

For j > 1 there exist three modes. One of them is 
purely magnetic, it has l = j, while two others are con¬ 
structed from the electric and longitudinal polarizations, 
each one of these two modes has an admixture of l = j + 1 
and l = j — 1 states. These two modes coexist for j > 1, 
while for j = 0 only one of them, which in this case has 
a purely longitudinal polarization and l = 1 is present. 

From Eq. one derives that the spectrum of the 
Coulomb problem is degenerate; it is triply degenerate 
provided n> j + 2 , j > 1 , doubly degenerate for levels 
with n = j + 1 , j > 1 , while the states which have either 
n = j or j = 0 remain non-degenerate. This conclusion 
agrees with the non-relativistic expansion, see Eq. ED- 

Interestingly, one and the same Sommerfeld formula 
Ea. (14.811 describes the discrete energy spectrum in the 
Coulomb problem for scalar, Dirac and vector particles. 
The only distinction is related to the angular momentum 
j, which takes the integer values j = 0 , 1 ... for bosons 
and half-integer j = 1 / 2 ,3/2 ... for fermions. 


V. CATASTROPHE WITH CHARGE 


Consider the charge density of a vector boson for a 
state with j = 0. Eqs. jttii . rrm give 


p w =2e 


t) 


e + ZZ- ) (v 2 + w 2 ) + - i (5-1) 


(5.2) 


where w defined by Eqs. irrra . iTTii . iiTTi reads 
1 


dv 2v 
£ + Zajr \dr r 
















In the region of small distances r <C Zajm Eo. 14.3911 
shows that ip oc r 7+1 / 2 . Consequently, from Eqs. 
E 22 ),E 1 we find the following estimates for v and w 


v 

w 


r 7 - 3/ 2 # 


1 + 1 / 2 r 7-3/2 

Za 


(5.3) 

(5.4) 


From here one derives an estimate for the charge density 
o> of the vector boson 


w 


— 2 e 


(l- 7 )(l + 27 ) 27 _ 4 
Za. 


r > 0 . (5.5) 


It diverges at the origin so badly that the total charge 
Q w = f p w d 3 r localized in any small sphere surround¬ 
ing the origin is infinite. 

The trouble does not stop here. Remember the density 
p = Z |e| S(r) of the Coulomb charge, which is located at 
the origin. This density results in the Y-term defined by 

En.iirrm 


t = 


ep_ 

m 2 


47 iZa 
m 2 


6{r) . 


(5.6) 


We did not consider it previously because the functions 
we dealt with were regular at the origin, allowing one 
to hope that their regular behavior makes the Y-term 
irrelevant. Since the charge density does not follow this 
pattern, we need to take the term given by Eq. 15.611 into 
account. The contribution of the (5-function in Eq. 15.61) 
to the boson charge density is given by the last term in 
Eci. 15.ID . which reads 


(p w ) T—term = w 2 ( 0) <5(r) . (5.7) 

m z 

Ea. 15.4D shows that w(0) = oo, which makes the density 
Eo. 15.7D infinite as well. 

We see that there are two closely located, though dif¬ 
ferent regions, which contribute to an infinite charge of 
the W-boson in the j = 0 state. One region is r > 0, 
where the density of charge Eq. (Oil behaves singularly 
as r —> 0. Another region is located strictly at the origin 
r = 0 , where an infinite coefficient w 2 ( 0 ) = oo in front 
of the (5-function in Eq. o makes the charge infinite as 
well. 

The origin of Eq. E3 can be traced down to the last 
term in Eq. I Ell . It contributes therefore to the charge 
density for all states. There is one more state, in which 
the coefficient in Eq. turns infinite, signaling a catas¬ 
trophic behavior of the charge. This is the “7 — 3/2” 
state with j = 1, see Eqs. (TOU.iTOTIi . To justify this 
statement, note that Eqs. (ra , im> and iTHTli imply 
that w oc r 7-3 / 2 . For j = 1 the inequality 7 < 3/2 
holds. Therefore for the state “7 — 3/2”, j = 1 one finds 
w( 0 ) = 00 , which makes the charge of the W-boson lo¬ 
cated strictly at the origin infinite. (There is no problem 
in that case with the charge in the region r > 0 .) 

The catastrophic behavior of the charge of the W- 
boson in j = 0 and j = 1 , “7 — 3 / 2 ” states was discovered 


in 5|, forcing the authors of this work to conclude that 
the pure Coulomb problem for IV-bosons is poorly de¬ 
fined. 


VI. VACUUM POLARIZATION 

Consider the conventional QED vacuum polarization. 
The potential energy of the W-boson propagating in the 
Coulomb field acquires an additional term, let us call it 
S(r ), which describes the polarization 

U(r) = -(l + S(r))— . (6.1) 

It suffices to consider the polarization effect in the lowest- 
order approximation, when it is is described by the known 
Uehling potential. Its small-distance asymptotic behav¬ 
ior is given by a simple logarithmic function, see e. g. 

m, 

S(r) ~ —a/3 In (mzr ), r —» 0 . ( 6 . 2 ) 

This function is related to the logarithm responsible for 
the scaling of the QED coupling constant 

a^ 1 ^) = a -1 (p 0 ) - f3 \n(p/p 0 ) . (6.3) 

The relation between Eqs. El and El is well-known, 
see e.g. book | 2 fl[, which presents it for one generation of 
leptons. The factor 0, which governs the scaling of the 
coupling constant and the potential in Eq. El equals the 
lowest coefficient of the Gell-Mann - Low /3-function. It 
is normalized here in such a way that for one generation 
of leptons 0 = (3 e = 2 / 37 r. 

It is important for us that a(p) rises with the mass 
parameter p, i.e. j3 is positive, /3 > 0; theoretical and 
experimental data agree on this fact, for a brief review 
see e. g. Ref. 24, the experimental data are provided by 
Refs. |23, 1211 |2j|. An estimation of /3 can be found from 
two reliable reference points a - 1 (m T ) = 133.498 ± 0.017 
and a~ 1 (mz) = 127.918 ± 0.018 provided in Ref. [2djl . 
Using them and taking the masses m T = 1776.99 + 0.29 — 
0.26 Mev and mz = 91.1876+ 0.0021 Gev recommended 
in J3j one derives from Eq. HO. 3D that 

0 ~ 1.42(1) . (6.4) 

More simple estimation of 0 can be done if one takes 
into account a contribution of all known charged fermions 
“naively” (neglecting complications, related to the QCD 
vacuum as well as possible contribution of scalars). This 
estimate yields 



Here summation runs over all charged fermions, qi is 
the charge of the fermion, the terms 1 and 3/5 in the 
bracket are due to the lepton and quark contribution for 
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one generation, the factor 3 after the bracket accounts 
for three generations. A discrepancy between “simple- 
minded” Eq. H6.5I) and more solid-based Eq. (IOI) is below 
20%. The normalization of the logarithmic function on 
the mass of the iT-boson rriz adopted in Eq. pre¬ 
sumes that the fine-structure constant a is taken at pre¬ 
cisely this scale, a = a(rnz) — 1/128. 

We are interested in high-momenta behavior in 
Eqs. iOl.iOl. where p, ~ 1/r m. An accuracy of 
Eas. lKi.4D . lKi.5D . as well as any other feasible estimation, is 
limited in this region by a contribution of unknown heavy 
charged fermions and scalars. However, this uncertainty 
does not affect our final conclusions. For our purposes it 
suffices to stick to a widely accepted hypothesis that (3 is 
a positive constant (or a slow-varying function up to the 
Grand Unification limit). 

Substituting Eqs. into Eq. one derives 


for r > 0, Ea. Ki.9D reduces to Ea. (l4.25D . Ea. Ki.9D can be 
rewritten in a more compact form 


v" + Gv' + Hv = 0, (6.10) 


where the coefficients G = G(r) and H = H(r) are 


G = 
H = 


2 U' | U' + V 
r + e-U + e-U - T ’ 

2 2 / U' U' + V \ 

r 2+ r\£-U + £-U- r) 


( 6 . 11 ) 

( 6 . 12 ) 


e-U-T 

e-U 


((e — U) 2 — m 2 ) . 


For a qualitative analyses it is convenient to eliminate 
the term with the first derivative by scaling the radial 
function v —> <p = <p(r) 


T(r) 


Za 2 /3 
m 2 r 3 


( 6 - 6 ) 


(£-U)(£- 


U - T) 



(6.13) 


where the lowest term of the a-expansion is retained. It 
is vital that for small distances, when r <C y/a/m, T(r) 
is positive and large, 


(This definition reduces to Ea. l4.26|) when T = 0). The 
classical equation of motion for W-bosons takes a simple 
form 


T(r) » \U(r)\ » m . (6.7) 

Note that the direct contribution of the vacuum polariza¬ 
tion given by the term S(r) in Ea. Ki.lD is not pronounced. 
In contrast, the T-term Eq. Ki.6D becomes dominant at 
small distances, making the effects related to the QED 
vacuum polarization very important. Since this term 
plays a crucial role below, let us verify its sign again. 
Consider a positive Coulomb center, Z > 0. Then the 
vacuum polarization produces negative charge density, 
p < 0. Since the charge of the W~ meson is negative, 
e < 0, we find from Eq. 12.19D that T = ep/m 2 > 0. We 
see that indeed, the T-term is positive, in accord with 
Ea. Ki.6D . 


A. Longitudinal polarization, j=0 

Ea, 14.24|l shows that a longitudinal state with j = 0 
is described by the single radial function v = v(r). 
Ea. I2.18D allows one to express the function w via v 

w = (e — U — T) _1 ( v' + 2 v/r) . (6.8) 

We need now to write the classical equation of motion for 
v, in which the term T is taken into account. The sim¬ 
plest way is to substitute W and w from Eqs. 
into Eq. RFM . which yields 

((e — U) 2 — m 2 )v = (6.9) 


-<p"+Uip = 0, (6.14) 

U = -H + G 2 /A + G'/2 , (6.15) 

where G , H are defined in Eqs. Ki. 1 lD . Ki.12D . Eo. Ki.14D 
can be looked at as a Schrodinger-type equation, in which 
U =U(r) plays the role of an effective potential energy. 

According to Eqs. K1.6D Ki.7D the T-term is large and 
positive at small distances. This fact makes the effective 
potential U (r) in Eq. K1.15D also large and positive when 
r —♦ 0 


U(r) ~ - H(r) ~ - U (r) T(r) ~ Z ® a ■ (6.16) 

Compare this result with the effective potential [U(r)]c 
for the pure Coulomb field. The latter one is a part of 
the Hamiltonian in Eq. (TP31) . For r —> 0 one finds from 
Eq. ll4.28D that 


[ ZT(r’) ] c — —(Za) 2 /r 2 . (6.17) 


It is taken into account here that the variables x and r 
in Eq. (TQ3I are proportional, see Eq. 14.27D . Eq. 16.16D 
shows that the vacuum polarization produces a strong 
repulsion in the effective potential M(r), in contrast with 
a mild attraction, which exhibits \U{r) ]c in Ea. Ki.17D for 
the pure Coulomb case. 

When the estimate Eq. is applicable, Eq. iKm 
allows an analytical solution 


ip(r) oc r exp 


Za (a/3) 1 / 2 


(6.18) 


It is taken into account here that Eq. ensures that 
V X W = 0. For a purely Coulomb case, when T = 0 


It shows that <p(r) exponentially decreases at small 
distances. According to Eqs. 16.8D . 16.13D the functions 
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v(r). w(r), also decrease exponentially here; correspond¬ 
ingly, the charge density of the W-boson Ea. H5.lll de¬ 
creases exponentially at the origin as well 


w 




a 1 

- 2 ex P 

m r z 


(a/3) 1 / 2 


Za (a/3) 1 / 2 


mr 


■ exp 


4 a 2 Zae 1 


m 


3T ex P 


Za (a/3) 1 / 2 


2Za (a/3) 1 / 2 


mr 


(6.19) 

( 6 . 20 ) 
. ( 6 . 21 ) 


Here a constant a depends on the normalization of v, 
which is specified in Eq. ra below. Eq. shows that 
for r > 0 in the vicinity of the origin the charge density 
is finite and small; which makes the charge located in 
this region finite as well. Eq. (IQ1 shows that u>(0) = 
0, which eradicates the contribution of the (5-function in 
Ea. (15.71) . Thus, the charge located strictly at origin r = 0 
is zero. 

We verified that an account of the QED vacuum po¬ 
larization erases the infinite charge of a vector boson for 
j = 0 state. 


B. Electro-longitudinal polarizations, j > 1 


where k satisfies 

k 2 = j 2 (j + l) 2 - 4(Za) 2 ^ 7 2 - 0 (j 2 - , (6.28) 

with 7 defined in Eq. 63- The two available values of 
A± should be attributed to the two electro-longitudinal 
modes. 

Comparing Eqs. (16.2511 . (16. 26D . which are valid when 
the vacuum polarization is taken into account, with 
Eos. iITTtI) and Eqs. emem which describe the 
purely Coulomb case, we see that the polarization 
changes drastically the behavior of the wave functions. 
One finds from Eqs. ffT27l . yr5BI that A± are positive for 
all j, j > 1, provided Z is not very large, Za < 1/2. 

From Eqs. 16.25ll . ()6.26D one deduces therefore that 
u(r),v(r ) —> 0, when r —> 0. Eq. El guarantees then 
that 'ic(O) = 0. As a result the contribution of the 15- 
function in Eq. ed to the charge density turns zero 
for all electro-longitudinal modes. This differs qualita¬ 
tively from the pure Coulomb case, which gives an infi¬ 
nite charge located at the origin for j = 1, = 0 state. 

We conclude that the QED vacuum polarization sup¬ 
presses the wave functions of a vector boson at the origin, 
eradicating thus the infinite charge of the boson, which 
plagues the problem for the pure Coulomb field. 


Eq.JO shows that in the region of small distances 
r <C a/m the T-term, which is related to the vacuum 
polarization, is large. This fact makes the function w in 
Eo. (12.1811 small, |w| <C \W\. As a result the asymptotic 
form of the equation of motion EH} at small distances 
reads 


w = V X (V X W) . (6.22) 

r z 

Ea. dOt ensures that VxlVis not zero identically pro¬ 
vided j > 1, which makes Eq. q> meaningful. 

Using Eq. (ET5l) to represent the electro-longitudinal 
modes and identities Eqs. o for the spherical vectors 
one rewrites Eq. i KTSU in terms of the radial functions 
u, v 


(Za) 2 u = —r 2 u" — tv! — u + y/j(j + 1) v , (6.23) 

(Za) 2 v = -y/j(j + 1 )(ru' +u- \Jj(j + 1) v). (6.24) 

Their solution is straightforward 

u = br x , (6.25) 

v = b yjj(j + 1) rX ' ( 6 - 26 ) 

Here 6 is a constant, and A can take one of the two pos¬ 
sible values, 


_ 1 j(j + 1 )±k 
~2 y 2 — 1/4 


(6.27) 


VII. NUMERICAL EXAMPLE 

To be more informative on the behavior of vector 
bosons in the Coulomb field let us solve the correspond¬ 
ing equations of motion numerically. Consider the j = 0 
state, describing it with the help of Eqs. (14.2411 and (IP1) . 
We need to specify the factor S(r), which describes the 
vacuum polarization in the potential in Ea. 16.ll l. For 
small r, r <C Za/m this factor plays a major role, 
while for larger it is less important. Let us construct 
a simple model, which gives a correct asymptotic be¬ 
havior Eq. dOi as r —> 0, and is physically reasonable, 
though not perfect, at larger r. Take with this purpose 
the Uehling potential, see e.g. 0, assuming that only 
charged leptons and quarks contribute to it 


Here 

F (x) = exp(— 2ir</) ^1 + ^ 1 dC > • ( 7 - 2 ) 

Summation in Eq. EH runs over all quarks and charged 
leptons, their charges qi and masses m„; are taken from 
Ref. |2J| . The model presented by Eq. O neglects com¬ 
plications related to the QCD vacuum, which may be 
substantial at large distances, but the role of the polariza¬ 
tion is insignificant in this region anyway. For small dis¬ 
tances r —> 0 the model Eq. ED reproduces the correct 
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xmc/fi rmc/A 




rmc !~h rmc/ti 


FIG. 1: The 2pi radial function v(r), which describes the 
n = 2, j = 0 discrete state of a VF-boson in the Coulomb 
field of a charge Z — 1. (a) Large distances r h/mc, v(r) 
is close to conventional non-relativistic wave function 2 p; (b) 
ultra-relativistic region r -C h/mc. Solid line - numerical 
solution, dashed line - analytical prediction of Eq. Klim . 


FIG. 2: The charge distribution p w (r) for the 2pi state 
of a IT-boson in the Coulomb field of a charge Z = 1. (a) 
Non-relativistic region of large distances r h/mc; (b) ultra- 
relativistic region r <§; h/mc. Solid line - numerical solution, 
dashed line - analytical prediction of Eq. em dashed-dotted 
line - the pure Coulomb case Eq. (TT31) . when the charge density 
diverges at the origin (shown with arbitrary normalization). 


asymptotic formula Eq. & in which the coefficient /3 is 
given by Eq. H6.5I) . A comparison of Eqs. dint and (liTfill 
shows that at small distances the accuracy of the model 
potential can be estimated as ~ 20%, which is sufficient 
for us. Using S'(r) from Eci. ( 17 . 11 ) one defines the poten¬ 
tial U(r) Ea. Ki.lD and the T-potential in Ecp ( 12 . 101 . After 
that the solution of Eq. (16.91) is straightforward. This so¬ 
lution should be normalized on the total charge e of the 
iF-boson 

e = Jp w (r)d 3 r, (7.3) 

where the charge density is defined by Eq. ED- 

Fig Cl shows the radial function v(r) for the 2pi state 
(n = 2, j = 0) for Z = 1. For large distances mr 1 the 
function v{r ) is close to the conventional non-relativistic 
wave function of the 2 p state, as it should be, see Section 
nm In the ultra-relativistic region mr <C 1 the func¬ 
tion v(r) changes sign, then shows an extremum, and de¬ 
creases exponentially when r —> 0, in accord with an ana¬ 
lytical estimate Eq. which is shown by the dashed 

line. The fact that the function v(r ) has a node in the 
relativistic region produces no controversy since v(r) is 


not a proper wave function and the conventional theo¬ 
rem, which counts the nodes of the wave functions for 
discrete levels is not applicable. 

In our discussion we did not try to construct proper 
wave functions, being content with a possibility to calcu¬ 
late the current. As an example, Fig[3 shows the charge 
density for the 2pi state in the Coulomb field of Z = 1. In 
the non-relativistic region mr 1 it behaves convention¬ 
ally. For the ultra-relativistic case mr -C 1 the density 
changes sign, exhibits an extremum and then decreases 
exponentially when r —> 0 in agreement with Eq. EHJ. 
Note that the “wrong” sign of the charge density, i. e. 
the positive charge density for the negatively charged W- 
boson (Fig[21 shows p w /e = — p H/ /|e|) produces no con¬ 
tradiction with general principles. The Pauli’s theorem, 
see e. g. Ref. [2(|, implies that the energy of a boson field 
is positively defined, but the sign of the charge density 
of the boson field remains not unequivocally determined 

m. 

The total positive charge located at small distances 
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proves be very small; for 2pi, Z = 1 state it is 

pro 

QW = 4tt / p w (r) r 2 dr ~ 0.676 ■ 10" 14 |e| . (7.4) 
Jo 

Here ro — 1.01 ■ lCU 3 m _1 is a node of p w (r). 

It is instructive to compare the found charge density 
with the one in the pure Coulomb problem, which is 
shown in Fig. [21(b) by the dashed-dotted line that repro¬ 
duces Eq. (1531) . One should be content in this case with 
an arbitrary normalization, since for the pure Coulomb 
field the normalization integral in Eq. (1731) is divergent. 
Fig. El (b) illustrates the fact that the vacuum polariza¬ 
tion reduces the charge density at the origin. 

The energy shift fe of the level 2pi (fe is a deviation of 
energy from the Sommerfeld formula Eq. iTTsli ) due to the 
vacuum polarization is found to be Se/m = —1.90 • 10 -7 . 
In relative units it is much bigger than the Lamb shift 
in atoms feLs/w e ~ Z 4 a 5 . The reason is obvious. The 
Lamb shift in atoms originates mostly from within the 
Compton distances r ~ r c = \/m e , which are smaller 
than the Bohr radius for the electron re ~ 1 /(Zam e ). 
For IT-bosons the situation is different. Light fermions, 
which contribute to Eq. o, allow the polarization po¬ 
tential to spread to large distances, as far as the Bohr 
radius of the W-boson, r ~ 1/rrq ~ 1 /(Zam). Therefore 
the energy shift due to the vacuum polarization gains 
substantial contribution from the non-relativistic region, 
where the wave function is large, which makes fe large as 
well (large compared to the Lamb shift in atoms). The 
accuracy of the energy shift calculations is limited by an 
accuracy of our model at large distances. The contri¬ 
bution of the QCD vacuum, which could be substantial 
here, is not described properly by the model based on 
Eq. (17.11) . Consequently, the presented above value for the 
energy shift should be considered only as an estimate. 

Nevertheless, one can derive an important lesson from 
this estimate. The found energy shift is small on the ab¬ 
solute scale, being lower than the non-relativistic binding 
energy by a factor of |&| x 4/ (Z 2 a 2 m) ~ 1.4-10 -2 . Thus, 
the dramatic variation of the function v{r) at the origin, 
which is produced by the vacuum polarization, makes 
only small impact on the spectrum. This is in contrast 
to a strong influence, which the vacuum polarization ex¬ 
ercises on the charge distribution of vector bosons. The 
fact that the energy shift is small makes the Sommer¬ 
feld formula Eq. EH a good approximation for discrete 
energy levels. 

VIII. DISCUSSION 

We demonstrated that the conventional QED vacuum 
polarization plays a very important, defining role in the 
Coulomb problem for vector bosons. Let us summarize 
the reasons leading to this conclusion. The Uehling po¬ 
tential, which describes the vacuum polarization in the 
simplest approximation is known to be a weakly attrac¬ 
tive and slowly varying function. For spinor particles it 


produces a small enhancement of the fermion wave func¬ 
tions on the Coulomb center. For vector bosons the sit¬ 
uation is different because the equations of motion for 
vector particles explicitly incorporate the external cur¬ 
rent. As a result, the density of the polarized charge 
p comes into the equations of motion for vector bosons. 
The corresponding term in the equations was called the 
T-potential, T = ep/rri 2 . The charge density p is nega¬ 
tive for an attractive Coulomb center, p < 0 when Z > 0, 
being singular on the Coulomb center, |p| ~ 1/r 3 . One 
derives from this that the vacuum polarization produces 
a repulsive T-potential, T = ep/m 2 oc 1/r 3 > 0 (remem¬ 
ber, e < 0). Since the T-terrn is singular at the origin, it 
plays a dominant role at small distances. 

Strong effective repulsion produced by the T-potential 
reduces the fields, which describe W-bosons on the 
Coulomb center. For j = 0 this reduction is dramatic, 
exponential. For j = 1, ”7 — 3/2” the suppression is of a 
more moderate power-type nature, but in both cases it is 
strong enough to eliminate the infinite charge, which is 
located at the origin in the pure Coulomb approximation. 

The above comments appeal to a chain of calculations. 
It is interesting to look at the obtained result from a more 
general perspective. The renormalizabity of the Stan¬ 
dard Model implies that by renormalizing relevant phys¬ 
ical quantities one is bound to obtain sensible physical 
results. The relevant quantity in question is the charge 
density of a vector boson. It follows from this that the 
important physical quantity, which should be renormal¬ 
ized, is the coupling constant. Its renormalization is ef¬ 
fectively fulfilled when the vacuum polarization is taken 
into account. Thus, it makes sense that the account of 
the vacuum polarization results in acceptable physical 
results. 

A proposed approach is very straightforward, which 
makes the Coulomb problem for vector bosons as simple 
and reliable as it is for scalars and spinors. All discrete 
energy levels can be easily evaluated, all relevant fields 
can be calculated and normalized properly. Presumably, 
all scattering data can also be evaluated, though the scat¬ 
tering problem was not discussed in detail in the present 
work. All these quantities include the Coulomb charge 
Z accurately, not relying on perturbation theory. Start¬ 
ing from this base, one can consider all other processes 
left outside the scope of the Coulomb problem by treat¬ 
ing them as perturbations. This includes the conven¬ 
tional QED processes, such as the radiative decay, pho¬ 
toionization, the radiative corrections. This includes also 
processes related to possible exchange of Higgs and Z- 
bosons. 

Previous attempts to formulate the Coulomb problem 
for vector bosons within the framework of the Standard 
Model have been facing a difficulty related to an infinite 
charge of the boson located near an attractive Coulomb 
center. This work finds that the polarization of the QED 
vacuum eradicates the problem. Usually the QED radia¬ 
tive corrections produce only small perturbations. It is 
interesting that in the case discussed the radiative cor- 
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rection plays a major, defining role. 
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APPENDIX A: HOMOGENEOUS MAGNETIC 
FIELD 

For a static homogeneous magnetic field Eq. (THU 
reads 

(e 2 - m 2 ) W = - (V - ieA) 2 W - 2 ieB x W. (Al) 

Assuming that the magnetic field is directed along the 
z-axis and introducing the new variables w a , a = 0, ±1, 
w± i = (W x ± iW y )/y/ 2, w 0 = W z one rewrites Eq. ED 
in a simple form 

(£ 2 — m 2 ) w a = — (V — ieA) 2 w a + 2eBaw a , (A2) 

which looks similar to the non-relativistic Schrodinger 
equation for a particle in the homogeneous magnetic 
field. This similarity allows one to write the spectrum 
Ea. (12. 2411 . 


APPENDIX B: RELATIVISTIC CORRECTIONS 
TO ENERGY LEVELS 


Here we present separate expectation values for four 
relativistic corrections in the same order as they appear 
in Ea. (13.611 . For l = 0 


j — 


i(Za) 4 



+ ° + - + ° 


(Bl) 


Here and below we specify the terms having zero expecta¬ 
tion values by writing the corresponding zeros explicitly. 
For l > 1 


4 E-nlj — 


i(Za) 


4 r 


rr 


[p 

1 ^ 

\8n 

21 + 1 J 


(Is) 


Z(Z + 1)(2Z + 1) 


+ 0 (B2) 


6(ls) 2 + 3(ls) -41(1 + 1) 

’Z(Z + l)(2Z-l)(2Z + l)(2Z + 3)_ 


where 


(ls) = ~\j(j + l)-l(l+l)-2] . 


(B3) 


Both Eo. dBlI) and Eo. (IB2l) lead to Eo. (l4.2ll . The to¬ 
tal relativistic correction would look very complicated 
and show no degeneracy if magnetic dipole or electric 
quadrupole moments of a vector particle are different 
from those values, which follow from the gauge theory. 


APPENDIX C: SPHERICAL VECTORS 


The conventional definition of spherical vectors, see 
e.g. Jfj, reads 

^ = v„iWvW+i), 

Yjrl = nY jm , (Cl) 

Yff=nxY$ . 

Here Yj m = Yj m (6,ip) is the spherical function, 
Y^,Yj£,Yj™' > are the electric, longitudinal and mag¬ 
netic vectors. The symbol V n in Eq. gd indicates the 
angular part of the gradient, VF(d,cj)) = V„F(0, <f>)/r, 
and n = r/r is a unit vector along the radius vector. 

Definitions Eqs. m imply the following properties of 
the spherical vectors 


• yS = -VW+ i) Yj 

(l) _ 
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V . y 1 ’" 1 — n 

(e) _ y(m) 
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V X Y% = -vW+l )4 
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1 jm — ^ 1 jm -> 
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V x Y K } — Y 

v A 1 jm — 1 jm 

r(l) 

" jm 


V x Y™ = -Y$ ^j(JTT) Yj 


(0 

jm 


(C2) 


(C3) 


The formulas for the Laplace operator read 

An Y$ = -j(j + 1) Y$ + 2vW+l) Y®, 

An Y$ = 2^/j(JTT) Y$ - (j(j + 1) + 2) 1 $, (C4) 
A n Y^=-j(g + l)Y^. 


Here A„ describes the angular part of the Laplacian, i.e. 
A F(d,(f>) = A n F/r 2 . The parity for electric and lon¬ 
gitudinal polarizations equals P = (-1)-', for magnetic 
polarization the parity is P = (—1) J+1 . The orbital mo¬ 
ment l takes the value l = j for the magnetic polarization, 
in agreement with the parity for this state. The electric 
and longitudinal polarizations are constructed as linear 
combinations of the two states with Z = j ± 1. For j = 0 
there exists only one spherical vector, which is purely 
longitudinal and has Z = 1. 


APPENDIX D: SPECTRUM OF 
ELECTRO-LONGITUDINAL MODES FOR j > 1 

Let us verify that for j > 1 the function g introduced in 
Eo. (14.2211 is nonzero. Consider first the “7— 1/2” mode. 
Substituting Eq. fl~T71 into Eq. IIT22I) one finds 

5 ~* a ^ (1/4 ~ 72)r7 ” 1/2 ’ (D1) 

which indicates that in this mode g is not zero. 

Consider now the ”7 — 3/2” mode, which incorpo¬ 
rates both possible polarizations at small distances. We 
need here the expressions for u and v at small dis¬ 
tances that are more accurate, then the ones in Eq. (TOI1 . 
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They can be derived by using mr <C 1 as a per¬ 
turbation in Eqs. 14.151) . 14. Kill . and pushing calculations 
one step beyond the simplest approximation given by 
Eos. 14.201) . 14.2111 . The result reads 


u —> b 1 — 


2 + (Za) 2 er \ ^_ 
7+1/2 Za) 


.7-3/2 


(D2) 


VJU+ 1 ) 


^7 — - — Zasr'j r 1 3 ^ 2 . (D3) 


Substituting Eqs. (Irnl.dmt into Eq. i TOH one finds that 
the main term oc r 7 ~ 3 / 2 cancels out in g , but the next 
one survives, giving 


9 ( 27 ~ 1 _ {Za)2 ^ r7 ” 1/2 ■ (D4) 


We verified that for j > 1 the function g is not an iden¬ 
tical zero for both electro-longitudinal modes. 


APPENDIX E: LONGITUDINAL MODE j = 0 
AND MATRIX MECHANICS 

In order to find the spectrum of the operator H in 
Ea. 14.2011 let us employ a method, which finds its in¬ 
spiration in an elegant approach to quantum mechanics 
developed by the Gotingen School and often called the 
matrix mechanics ; the book Ref. [23| gives its systematic 
presentation. We modify it for our purposes as follows. 
Assume that one needs to find discrete spectrum of some 
Hermitian operator 7 H. (in our case it is the operator Ti 
in Eq. (IP3l) h Let us presume that one is able to find the 
operator 9 , which satisfies 

n = e^e + \ Q , (ei) 

where Ao is a number. Define then a new operator TL, 

n = ee^ + \ 0 . (E2) 


Let us verify that the two operators 7+ TL have very sim¬ 
ilar sets of eigenvalues. Consider an eigenfunction p of 
TL , with the eigenvalue A 


TLp = A ip . 

(E3) 

Taking 


P = 9p , 

(E4) 

one verifies that 


TL <p = {9 fit + A 0 )9 p = 9(9U + \ 0 )<p 
= 9 TLp = \ 9 p = Xp . 

(E5) 

This shows that either the function p is an eigenfunc¬ 
tion of TL with the eigenvalue A, or p = 0. The first 
options makes A an eigenvalue of both operators TL, TL. 
The second one implies that 


(E6) 


which indicates that A = Ao is a candidate for an eigen¬ 
value of TL because Eq. implies TL p = XqP- Eq. dt 
provides a convenient way to derive the corresponding 
eigenfunction. There is though a subtlety here. The 
found from Eq. 1E6I) p may, or may not satisfy the bound¬ 
ary conditions. If it does, then it represents the eigen¬ 
function and A = Ao is an eigenvalue. Otherwise, Ao does 
not belong to the discrete spectrum, as would be the case 
in an example discussed. One should also verify that an 
action of the operator 9 in Eq. dl> (as well as the op¬ 
erator in Eq. 1E8I) below) does not spoil the boundary 
conditions. We presume here that this is the case, and 
verify later on that this assumption holds for a particular 
example discussed, see Eq. 

We conclude that any discrete eigenvalue of TL is an 
eigenvalue of TL as well, with one possible exception of 
Ao- By reversing the argument, one derives that if p is 
an eigenfunction of TL with the eigenvalue A, 


TLp = A p 


(E7) 


then 


ip = 9^ip 


(E8) 


satisfies Eci. (IE81) with A = A. We see that the two sets of 
discrete eigenvalues of the two operators TL,TL are same, 
except for Ao, which may, or may not be present in one, 
or both sets of spectra. The crucial for us point is that 
the operator TL can be more simple for analyses than the 
initial operator TL. 

Taking the operator TL from Eq. , we construct 
the operators 9 , 9 1 in the form 


9 = 
9 f = 


d_ 

dx 

d 

dx 


b c 

— 1 -T7 ’ 

x x + 1 

b c 

- 1 - 77 > 

x x + 1 


(E9) 

(E10) 


where a,b,c are real numbers. From Eas. (IE9|) . IE10I) it 
follows that 


e<e = -£+J + 2b(a + c) + 6(6- 1) 

dx z x xT 

2 c(a — b) c(c — 1) 

x + 1 (x + l) 2 

There are four x-dependent rational terms in Eq. cm 
while only three coefficients a, b , c are available for tun¬ 
ing to make them identical to similar terms present in 
the operator TL. However, the coefficients in Eq. iron 
prove to be “user-friendly”, making this procedure pos¬ 
sible. Taking 


a = b = 7+1/2, c = —1 , (E12) 

Aq = -(7 + 1/2) 2 , (E13) 

one satisfies Eq. & Taking 9 , 9 1 defined in 
Eas. IE9I) . IE10I) and (IE 121) one constructs TL Eq. llE5ll . 


9 ip = 0 
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with the result given in Eq. g31. The “nasty” singular 
at x = — 1 term disappears from Ti. The latter opera¬ 
tor describes a conventional Coulomb-type problem with 
L = 7 + 1/2 playing a role of an effective (non-integer) an¬ 
gular momentum. From Eq. <TX~32l) one finds that regular 
at x = 0 solution of the eigenvalue problem Tip = —x 2 p , 
satisfies Eq. 14.341) . Ea. 14.3611 . which ensures that this so¬ 
lution is regular at infinity, completely defines a set of 
discrete eigenvalues of H. 

The set of eigenvalues of Tt gives the eigenvalues of 
the original operator 7 i, except for possibly one addi¬ 
tional eigenvalue Ao, which is discussed below. The 
eigenfunctions of Ti can be found from Eq. lESIl . Us¬ 
ing Eqs. iTeTTI . IET 2 I 1 one presents them in a form of 
Ea. 14.3111 . 

In order to verify whether Ao is an eigenvalue of Ti one 
needs to find + from Ea. 1E6I) . Ea. 1E9I) gives 

(~~J~ + (7 + !/ 2 ) ——— r') V 5 = 0 , (E14) 

\ ax x x +1 / 


which leads to 

p = (x + 1 ) _1 x 1+1/2 exp [ (7 + 1/2) a;] . (E15) 

Since this function is singular at x = 00 , it cannot be an 
eigenfunction. Consequently Ao is not an eigenvalue. 

The function ip defined by Eq. (TPT1) exhibits regular 
behavior at both boundaries Eqs. (ipgii . trogi) . This en¬ 
sures that p is an eigenfunction. Note, that specifying the 
operators 9, 9 f one had an additional option. One could 
have chosen in Ens. lE 1 21) and all the following relevant 
formulas —7 instead of 7 . It this case, however, instead 
of Eci. 14.31)11 one obtains p oc a ; -7+1 / 2 for x —> 0, which 
indicates a singular, unacceptable for an eigenfunction 
behavior. 

We conclude that the full set of all discrete eigenval¬ 
ues of Ti is specified by Eq. ( 10 ( 1 . The corresponding 
eigenfunctions are given by Eqs. (10111 . (1071 ). 
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